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LONG-TERM  GOALS 

My  long-term  goal  is  improving  our  understanding  of  interactions  between  sub-grid  scale  models  and 
numerical  errors,  which  affect  simulation  skill  for  physical  and  biological  processes  in  estuaries  and 
coastal  ocean  environments.  For  prognostic  simulations,  slow  accumulation  of  truncation-level  errors 
will  progressively  distort  model  results.  Present  sub-grid  scale  models  selectively  address  certain 
processes  (e.g.  turbulence)  while  neglecting  other  processes  (e.g.  small-scale  internal  waves, 
particulate  nature  of  algae).  As  small-scale  nonlinear  errors  can  accumulate  over  time,  understanding 
the  errors  and  modeling  the  neglected  processes  are  critical  issues  in  improving  our  prognostic  skill. 

OBJECTIVES 

I  am  examining  the  un-physical  behavior  of  internal  waves  modeled  under  the  hydrostatic 
approximation.  This  knowledge  is  a  precursor  to  developing  improved  modeling  techniques  that 
account  for  transfer  of  energy  from  resolved  internal  waves  into  sub-grid  scale  internal  waves.  These 
higher  frequency  waves  can  break  on  sloping  boundaries  or  lead  to  enhanced  wave-wave  interaction, 
thereby  amplifying  vertical  mixing  and  changing  the  spatial  and  temporal  evolution  of  the  oceanic 
density  structure. 

APPROACH 

Nonlinear  terms  in  the  momentum  equations  slowly  cause  an  initially-linear  internal  wave  to  steepen, 
resulting  in  nonlinear  and,  what  should  be,  non-hydrostatic  evolution.  However,  in  a  hydrostatic 
model,  the  non-hydrostatic  dispersion  is  missing,  therefore  modeled  internal  waves  steepen  until 
numerical  diffusion,  dispersion,  or  dissipation  provides  balance.  Thus,  steepening  internal  waves  in  a 
hydrostatic  model  will  artificially  diffuse  the  density  gradients,  disperse  the  waves  or  dissipate  their 
energy.  Alternatively,  a  theoretically  error-free  hydrostatic  model  must  result  in  waves  that  always 
steepen  until  breaking  -  thus  artificially  diffusing  the  density  gradient  in  breaking-induced  mixing. 

We  are  examining  approaches  for  modeling  the  correct  steepening/dispersion  balance  without  resorting 
to  the  computational  expense  of  a  full  non-hydrostatic  solution.  I  am  working  with  two  University  of 
Texas  (UT)  graduate  students,  Ms.  B.  Wadzuk  and  Ms.  S.  Delevan,  applying  hydrostatic  and  non¬ 
hydrostatic  models  of  the  Navier-Stokes  equations  to  simulate  internal  waves  in  simple  geometries. 
These  results  will  be  compared  to  laboratory  experiments  and  the  KdV  wave  model  of  Horn  et  al  2000. 
After  validating  the  models,  we  will  examine  the  divergence  of  results  between  the  hydrostatic  and 
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non-hydrostatic  models  under  varying  model  eonditions  (e.g.  grid  resolution,  time  step,  stratifieation). 
By  separately  quantifying  errors  in  diffusion  of  mass  and  momentum,  we  will  look  for  a  means  of 
predieting  required  model  parameters  for  adequately  resolving  internal  waves  with  a  hydrostatie  eode. 
In  the  past  year,  work  on  separating  and  analyzing  numerieal  diffusion  has  been  done  in  eollaboration 
with  Asst.  Prof.  B.  Laval,  University  of  British  Columbia  (UBC)  and  Prof  J.  Imberger,  University  of 
Western  Australia  (UWA).  Collaboration  with  Prof  Imberger  also  involves  development  of  nested 
grid  methods,  whieh  are  neeessary  for  higher  resolution  where  internal  waves  are  breaking.  One  of  my 
UT  graduate  students,  Mr.  J.  Fumans,  is  presently  at  UWA  on  a  Fulbright  Seholarship  and  is 
eondueting  this  work. 

WORK  COMPLETED 

This  first  full  year  of  the  projeet  has  been  foeused  on  eompleting  preliminary  work  for  quantifying 
numerieal  diffusion,  while  simultaneously  building  student  expertise  and  involvement.  This  projeet’s 
first  graduate  student,  Ms.  Wadzuk,  has  eompleted  her  M.S.  degree  and  has  begun  Ph.D.  studies.  Her 
M.S  researeh  report,  entitled  Evolution  of  Internal  Waves  in  Hydrostatic  Models:  A  Study  of  Dynamic 
Pressure,  is  available  as  a  web  download  at  the  projeet  web  site  listed  above.  Ms.  Wadzuk  is  in  the 
midst  of  developing  a  praetieal  approaeh  for  implementing  the  solution  of  non-hydrostatie  pressure  in 
the  framework  of  a  hydrostatie  numerieal  model.  We  have  eompleted  the  literature  review  and 
analysis  of  prior  approaehes.  In  the  past  year,  Ms.  Wadzuk  has  developed  profieieney  in  Fortran  95 
eoding  and  investigated  Poisson  methods  for  the  solution  of  a  fraetional-step  equation  as  part  of  the 
non-hydrostatie  solution.  A  seeond  graduate  student,  Ms.  Delavan,  joined  this  projeet  over  the  past 
summer.  Her  foeus  is  on  eondueting  and  analyzing  simulations  to  quantify  the  performanee  of  the 
hydrostatie  model  over  a  range  of  internal  wave  eharaeteristies.  During  the  summer,  Ms.  Delavan  has 
developed  the  non-dimensional  formulation  of  the  numerieal  experiments  and  beeome  familiar  with 
applying  the  hydrostatie  model. 

With  Asst.  Prof.  Laval  (UBC)  and  Prof.  Imberger  (UWA),  we  have  eompleted  a  study  on  quantifying 
and  redueing  the  effeets  of  numerieal  diffusion  in  models  of  internal  waves  (Laval  et  al.  2002).  In 
related  eollaborative  work  with  Dr.  Chris  Dallimore  (also  of  the  UWA),  we  have  shown  that  2D  and 
3D  models  ean  be  linked  to  improve  the  representation  of  the  propagation  of  density  eurrents  down  a 
slope,  whieh  has  resulted  in  two  more  papers  aeeepted  to  the  Journal  of  Hydraulic  Engineering  this 
past  year.  This  linking  of  models  of  different  spatial  dimensions  provides  another  avenue  for 
investigating  improvements  to  internal  wave  modeling.  We  have  developed  a  hypothesis  that  an 
inviseid  wave  model  for  a  2D  density  surfaee  eould  be  linked  to  a  3D  hydrodynamie  model  as  a  means 
of  improving  the  internal  wave  representation. 

RESULTS 

The  first  diffieulty  in  examining  the  evolution  of  internal  waves  in  numerieal  models  is  separating  the 
effeets  of  numerieal  diffusion  of  mass,  numerieal  dissipation  of  energy,  and  numerieal  dispersion  of 
waves.  We  have  learned  that  we  ean  quantify  and  eompensate  for  numerieal  diffusion  by  eomputing 
ehanges  in  baekground  potential  energy  oeeurring  during  adveetion.  This  builds  on  eoneepts  from 
prior  ONR- funded  studies  of  Winters  et  al.  (1995).  It  has  beeome  elear  that  quantifying  and 
eontrolling  numerieal  diffusion  is  a  neeessary  preeursor  to  evaluating  the  steepening  of  internal  waves 
under  the  hydrostatie  approximation.  Unless  the  numerieal  diffusion  of  mass  is  limited,  the  steepening 
of  waves  will  be  eontrolled  by  numerieal  diffusion  rather  than  by  the  governing  equations.  The  result 
is  degradation  of  the  thermoeline  gradients  and  distortion  of  wave  evolution.  In  a  paper  reeently 
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accepted  by  the  Journal  of  Hydraulic  Engineering,  we  demonstrate  an  approach  applying  a  sharpening 
filter  to  the  density  field.  This  ereates  an  a  posteriori  adjustment  that  eompensates  for  numerical 
diffusion’s  effect  on  the  density  strueture.  This  has  been  demonstrated  (Figure  1)  to  provide  a  better 
representation  of  the  evolution  of  the  sharp  density  structure  as  seen  in  field  data.  In  other  words, 
standard  application  of  a  numerical  model  leads  to  slow  numerical  diffusion  of  the  thermoeline,  whieh 
distorts  the  internal  wave  evolution.  By  quantifying  and  compensating  for  this  diffusion  with 
sharpening  filter,  we  can  retain  a  better  model  of  the  thermoeline  gradient  and  improve  the  model  skill. 
The  major  difficulty  with  the  approach  is  that  it  is  inherently  ad  hoc  -  i.e.  it  is  an  artificial  method  of 
eountering  error  after  it  has  occurred  rather  than  preventing  it  during  advection.  The  method  eannot 
predict  precisely  where  numerieal  diffusion  has  oeeurred,  but  only  the  net  amount  over  the  whole 
domain,  whieh  then  is  removed  by  artificially  sharpening  the  thermoeline  throughout  the  system. 

Thus,  this  has  some  parallels  to  “global  fill”  approaches  that  are  used  to  globally  redistribute  the  net 
mass  error  in  models  with  non-conservative  transport  routines.  However,  instead  of  being  a  global 
redistribution  of  mass,  this  is  a  global  readjustment  of  density  gradients  (with  weighting  at  the 
thermoeline).  This  filtering  method  has  been  sueeessful  in  controlling  increases  in  numerieal  diffusion 
over  the  course  of  a  simulation  (Figure  2). 


Figure  1.  Temperature  gradient  v.  depth  for  field  data  and  model  results  with  and  without  filter, 
[graph:  Field  data  temperature  gradient  and  filtered  model  results  show  peak  at  -2.5  C°cm'^  while 

unfiltered  model  results  has  peak  at  -1.2  C°cm'\j 
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Figure  2.  The  change  in  background  potential  energy  (Eb)  for  filtered  and  unfiltered  model  results 
(normalized  by  the  change  in  the  filtered  result).  These  are  simulations  of  adiabatic,  periodically 
wind-forced,  stratified  basins  for  which  a  perfect  hydrostatic  model  would  show  zero  change  in  Eb. 
[graph:  The  unfiltered  Eb  results  show  an  increase  15  times  greater  than  the  increase  in  the  filtered 

results  over  the  course  of  a  year.] 
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IMPACT/APPLICATIONS 


Improvements  in  our  ability  to  eorrectly  represent  internal  waves  in  hydrostatie  models  will  provide 
improved  representation  of  the  spatial  and  temporal  evolution  of  the  density  field  in  the  oeeans.  This 
will  improve  our  ability  to  use  model  results  to  prediet  thermoeline  position,  thiekness  and  density 
gradients  -  espeeially  in  eoastal  regions  where  steepening  and  breaking  of  internal  waves  along  the 
slope  may  be  an  important  source  of  mixing.  The  approach  developed  in  Laval  et  al.  (2002)  provides  a 
practical  method  of  evaluating  the  magnitude  of  numerical  diffusion  and  reducing  its  effects. 

TRANSITIONS 

The  background  potential  energy  algorithm  and  density-filtering  algorithm  of  Laval  et  al.  (2002)  has 
been  implemented  in  distribution  version  1.5.2  of  the  ELCOM  model  (Hodges,  et  al.  2000).  In 
addition  to  the  present  research  program,  the  model  is  being  used  at  Stanford  University,  University  of 
Stuttgart,  University  of  Western  Australia,  University  of  British  Columbia,  Instituto  Nacional  Del 
Agua  Y  Del  Ambiente  (Argentina),  Kinneret  Limnological  Laboratory  (Israel),  University  of  loannina 
(Greece). 

RELATED  PROJECTS 

The  long-term  improvement  of  models  for  sub-grid  scale  effects  is  the  subject  of  ongoing  collaborative 
efforts  of  the  PI  with  Prof.  Jorg  Imberger  (Stockholm  Water  Prize  Laureate,  1997)  at  University  of 
Western  Australia.  These  efforts  include  development  of  a  benthic  boundary  layer  underflow  model, 
development  of  particle  tracking  algorithm,  and  implementation  of  nested  grids  for  increased  local 
model  resolution.  A  University  of  Texas  graduate  student,  Mr.  Jordan  Fumans,  who  is  under  my 
supervision  but  is  presently  studying  at  the  University  of  Western  Australia  with  a  Fulbright 
scholarship,  is  accomplishing  the  work  on  nested  grids.  This  work  is  directly  related  to  the  study  of 
internal  wave  evolution,  as  we  would  like  to  apply  higher  grid  resolution  in  regions  where  waves  are 
rapidly  evolving  or  breaking. 
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